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Abstract 

We show that the existence of a possible fourth heavy generation of quarks gives 
rise to a significant enhancement to the neutron electric dipole moment in comparison 
with the Standard Model prediction. The smaller degree of suppression in this case 
is linked to the presence of the operators of dimension < 6 which enter into the 
effective Lagrangian with coefficients proportional to the square of the top quark 
mass. Numerically, the enhancement is mainly associated with chromoelectric dipole 
moment of the s quark which appears at three loop level, of the order a s a^m s m|/m4 
from the CP-odd combination of mixing angles between second, third and fourth 
generations. Its value is calculated explicitly in the limit of large masses of the fourth 
generation of quarks. The corresponding contribution to the electric dipole moment 
of the neutron is 5 • 10 _30 e • cm in the most optimistic scenarios about the values of 
the Kobayashi-Maskawa matrix elements. The additive renormalization of #-term in 
this model is estimated as 1CP 13 . 
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1 Introduction 



The Kobayashi-Maskawa (KM) model looks now as the most natural description of CP- 
violation. It describes properly CP-odd phenomena in the decays of neutral i^-mesons and 
predicts extremely tiny CP-odd effects in the flavour- conserving processes. The current 
experimental limit on the electric dipole moment of neutron (EDM) 

d N /e < 1(T 25 cm, (1) 

exceeds the realistic Standard Model prediction for this quantity by seven orders of magni- 
tude. However, this gap between theory and experiment is of great use for limiting a new 
CP-violating physics beyond the Standard Model. The purpose of this work is to consider 
the electric dipole moment of neutron in the model with an additional heavy generation of 
quarks preserving the same KM origin of CP-violation. 

The reason for introducing an extra heavy generation of quarks into physics comes from 
a relatively large mixing in B — B meson system 0. Its existence is, of course, questionable. 
However, the set of constraints on mixing angles and unknown heavy masses can be derived 
from the low energy phenomenological data. The analysis shows that the existence of 
additional heavy quarks with masses not lighter than m t is not excluded. We would use the 
conclusions developed in the work for possible values of KM matrix elements in terms 
of the Wolfenstein parameter A = \V US \ = 0.22. 

Let us denote the fourth generation flavours as (h g). Then the best scenario about a 
large mixing between the third and fourth generations consistent with current experimental 
data is given by: 

\V hb \ ~ \V tg \ ~ 0(A) (2) 

This scenario is quite natural if we assume the masses of h and g quarks lying below the 
perturbative unitarity limit of 500 GeV. Using the unitarity conditions for the KM matrix 
and already known values of matrix elements we deduce the following: 

\Vcg\ ~ |T4 S | ~ 0{\ 2 )- \V hd \ ~ \V ug \ ~ 0(X 3 ). (3) 

The enlarged KM matrix of this model possesses 9 independent parameters which 
include six mixing angles and three CP-violating phases. To avoid the uncertainties of 
reparametrization, we would describe all CP-odd flavour-diagonal amplitudes in terms of 
imaginary part of three independent quartic combinations of KM matrix elements: 

\m{V: d V tb V; b V cd ) ~ 0(X 7 ); Im^V^V^) ~ 0(X 7 ); 

lm(V t :V tb V: b V hs ) ~ G(A 5 ) (4) 

All other rephasing invariants could be reexpressed using these three combinations and 
moduli of KM matrix elements only. 
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2 Standard Model prediction 



Let us start from the Standard Model (SM) prediction for the electric dipole moment 
(EDM) of neutron. The violation of the CP-symmetry in the Standard Model originates 
from the sole complex phase in the KM matrix. To lowest, quadratic order in the weak 
interaction all CP-odd flavour-conserving amplitudes turn to zero trivially. The point is 
that in this approximation those amplitudes depend only on the moduli squared of elements 
of the KM matrix, so the result cannot contain the CP-violating phase. In the next, quartic 
order in semi- weak coupling constant g w , the expression for EDMs of quarks vanishes at 
two- loop approximation without hard gluon radiative correction taken into account ||. It 
can be shown, however, that the inclusion of one hard gluon loop prevents EDMs from 
identical cancellation. 

The values of EDMs to one gluon loop accuracy was calculated first by Khriplovich Q]. 
We quote here his result: 

d * = e6 ina 5 lo S —2 lo S -| - e • 2 • 10- 34 cm, 5 
1087T J mf m l c 

where Gp = y/2g^/(8M 2 ) is the Fermi constant and is the mass of i-flavoured quark. 
S denotes here the only possible CP-odd invariant of 3 by 3 KM matrix. The analogous 
result for d u is much smaller being proportional to m 2 s . 

Formula ([]) was obtained in the limit of the effective four-fermion contact interaction 
which is valid when all quark masses and characteristic loop momenta are much smaller 
than masses of SU(2) gauge bosons. This is not exactly true in general because the top- 
quark does not satisfy this demand. However, we could use the expression (||) as a good 
estimate for EDM of d-quark in the Standard Model. The reason being is that m t enters 
to (|5p under logarithm only and may be replaced by m w . 

The EDM of neutron in Standard Model, however, is much larger than EDMs of its 
constituents due to the so called "long distance" effects [§, §, [7]. The most reliable estimates 
according the rules of chiral perturbation theory predicts the EDM of neutron at the level 
djv — 2 • 10 _32 e ■ cm. The two orders of magnitude enhancement here is basically due to a 
smaller number of closed loops and bigger factor connected with strong interactions. 

Our basic idea is to extend this consideration to the case of KM model with an additional 
heavy generation of fermions. The reason is that instead of one CP-violating phase, the 
enlarged variant of this model possesses three CP-odd parameters and additional flavour's 
combinations in the amplitudes. It may cause the contributions to EDMs being propor- 
tional to the square of the mass of the top-quark. As a result, short distance contributions 
to EDMs of quarks would be suppressed not by G 2 F m 2 but rather by a w Gp. 

In general, the problem of EDM calculation should be divided into two independent 
parts. First, we construct a low energy effective Lagrangian in terms of u, d, s quarks, 
gluons and external electromagnetic field. Then, we recalculate this Lagrangian to the 
EDM of neutron using all available methods for doing low energy hadronic physics. The 
second part is independent from the concrete model of CP-violation at high energies. It 
is clear that in our case the possible enhancement should be associated with effective 
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operators of low dimension, not bigger than 6. It is easy to identify all these possible 
structures. Integrating out c, b, t, g, h quarks and SU(2) gauge fields, we construct a low 
energy effective Lagrangian in the form: 

i i 

(6) 

where (Fa) denotes F^j^ and the summation is held over light flavours: i = u, d, s. The 
first term in (|6]) represents the induced #-term, perturbative contribution to the total #-term 
of the theory. The next operator of dimension 6 was introduced originally by Weinberg 
||; the SM value of corresponding coefficient Cw was calculated in ||. Other terms in 
(§) with dim=5 are the operators of EDMs and chromoelectric dipole moments (CEDM). 
The dimension 5 of these operators in some sense is fictitious because of a chirality flip 
making coefficients di to be proportional to light masses m^. Therefore, both di and di are 
suppressed by at least two powers of heavy mass corresponding to the weak interaction scale. 
Other CP-odd operators of dimensions higher than 6 are unimportant in our consideration 
because they are suppressed by additional powers of heavy masses. 

Some comments should be made at this point. We do not add to (0) CP-odd mass 
operators of quarks iq^qi because they could be incorporated to the #-term by mean of 
the chiral rotation. Other flavour conserving CP-odd operators of dim=6 built from four 
quark fields are suppressed in comparison with CEDMs. The last fact refers to our concrete 
model of CP-violation. Because of the V - A character of the theory, effective operators are 
originally formulated in terms of left handed fields. In its turn, CP- violation may only arise 
through the chirality flip which gives extra powers of light mass and affects the suppression 
of four quark operators. This last remark is connected with the so called "axial polyp 



operators" introduced originally in JLO]: 



cpqan{G^,D u }j 5 qi. (7) 

This operator is identical to CEDM one if we assume the equations of motion to be pre- 
served: 

ij^D^qi = miqi (8) 

Further analysis is recalled to single out leading operators in (||) and find corresponding 
coefficients. 



3 Flavour's structure of C e ff 

We start investigating the coefficients in (^) by determining the flavour's arrangement along 
a fermion line. Let us denote by f the Green function of f-flavoured fermion. Then a CP-odd 
amplitude for fermionic operators in quartic order in semi-weak constant could be written 
in the following form: 

E iMv; f v Jk v; k v lf ) mf. (9) 

3,k,l 
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For pure gluonic operators the corresponding structure looks as: 

E Hva{V* f V jk V: k V lf ) fjkl, (10) 

fj,k,l 

where the cyclic permutation of the kind fjkl = Ifjk = klfj = jklf is allowed. 

In three generation formulation of KM model one has the only source of CP violation 
which fixes unique flavour structures for both type of operators [[§]. In the four family 
case, the number of CP violating phases is three. Their concrete choice could be done in 
a different ways complicating our analysis. However, it does not affect the main property 
responsible for the cancellation of EDMs at two loops. It is easy to see that independently 
on the number of families the expression @ is antisymmetric under the interchange of 
flavours j and 1: 

£ i\m{V; f V jk V lk V lf ) fjklf = i £ ihn{V* f V jk VC k V lf ) f(jkl - lkj)f. (11) 

j,k,l j,k,l 

This antisymmetry is sufficient to set to zero both EDM and CEDM of quark to two-loop 
approximation |H|. 

It is useful to classify a variety of CP-odd amplitudes by mean of dynamical arguments. 
The enhancement of the short distance contributions to the £ e ff which we expect to get 
is closely related to the fact that all characteristic loop momenta are comparable with 
the weak interaction scale. Therefore, inside the loops, we are legitimate to put all quark 
masses to zero except m t , m g and m^. In other words, inside the loops, we are able to 
identify propagators of light quarks: 

c = u = U] d = s = b = D. (12) 

It should be mentioned here that c and b quarks play a twofold role. They are considered 
as heavy quarks at normal hadronic scale and almost massless inside loops, at the scale 
of weak interactions. This requirement immediately leads to the vanishing of pure gluonic 
operators as well as EDM and CEDM of u-quark. Indeed, applying the unitarity condition 
for the KM matrix 

VlV dk + V; u V sk + V b * u V bk = 6 uk - V* u V gk (13) 
we perform the summation over j and / in (|9|) explicitly: 

£ Hm(V^V jk V lk V lu ) ujklu = ]T ilm( V gu V gk V gk V gu )u(Dkg - gkD)u = 0. (14) 

j,k,l k 

For s and d quark operators the situation is quite different. An analogous procedure 
for them leads to the following combination: 

- tim{ v t*f V tkVh k V hf )f[tkh - hkt + Ukt - tkU + hkU - Ukh}f, (15) 
2 k 

where / = d, s. 

Now we will concentrate on s-quark operators because its mixing with third and fourth 
generations is a priori bigger than that of d-quark. 
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Taking the last sum over k we obtain: 

(16) 

The rephasing invariant combination of KM matrix elements in (|T6| ) to a good accuracy 
coincides with that responsible for CP-odd Bg meson mixing: 

Im(V*V tg V h * g V hs ) = -lm(V*V tb V: b V hs ), (17) 

and the resulting expression takes the form: 

l -hn{V:y th V^ b V hs )s[t{b-g)h-h{b^ (18) 

It may be of the order A 5 and it makes the electric and chromoelectric dipole moment of 
s-quark the most important between other CP-odd operators. 

For the d-quark all considerations presented above are valid with the only replacement 
in rephasing invariants: 

lm(V t :V tb V^V hs ) — > lm(V; d V tb V: b V hd ) ~ 0(X 7 ). (19) 



4 EDM of neutron 



When the flavour structure is fixed, it is possible to find the relative meaning of different 
operators for the EDM of neutron using order of magnitude estimations for corresponding 
coefficients. 

First we take an electric dipole moment of d-quark. Combining together all phase space 
factors, coupling constants and taking into account (fl9|) we get the following estimate: 

d d ~ A 7 ^-^m d 4- (2°) 
47T Air lbiT z m A w 

In the next section, we will find that this estimate is reasonable. Now using a simplest 
constituent model, we obtain the corresponding contribution to the EDM of neutron at the 
level: 

d N ~ d d ~ e • 3 • 10~ 32 cm. (21) 

It turns out that it does not exceed the Standard Model prediction. The enhancement in 
comparison with SM value of d d is just two orders of magnitude instead of m\fm\ ~ 10 4 . 
It can be explained rather trivially. In contrast with the SM result (|^), the estimate (PH|) 
does not possess any logarithmic enhancement. Another origin of this deficiency is related 
to a smaller value of a s in ( PU| ) and smaller numerical factor in comparison with that of 



In the same manner we estimate the chromoelectric dipole moment of s-quark: 

,a s a w 1 m\ 5 a s a w 1 m\ 

d s ~ lm{V ts V tb V hb V hs ) — — ——G F m s — ~ A ^^Trr^G F m s — • 22 

47T 47r lbTC 2 mi, 47r Air Ibn 1 m l 
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The renormalization to hadronic scale does not seriously change this estimate. The recal- 
culation of the contribution of this quantity to the EDM of neutron is a separate problem. 
Here we use the result of Khatsimovsky, Khriplovich and Zhitnitsky |L1] which does not 



exhibit any additional suppression of the contribution to NEDM from this operator: 

d N ~ --d a ~ e • 5 • 10" 30 cm. (23) 

This estimate shows that the four-generation formulation of the KM model in the most 
optimistic scenario about possible values of CP-odd rephasing invariants leads to the electric 
dipole moment of neutron two orders of magnitude bigger than its SM value. The electric 
dipole moment of s-quark appears to be of the same order of magnitude as d s . However, its 
contribution to NEDM seems to be suppressed in comparison with that of CEDM operator 
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It is useful also to estimate values of other terms in £ e ff- Both the induced #-term and 
Weinberg operator turn out to be suppressed by the ratio m\jm 2 w . It means also that the 
characteristic loop momenta could range widely, from m& to electroweak scale. Therefore, 
it is quite possible that these contributions would match also a large logarithmic factor so 
the total suppression would be of order m\jm 2 w log (rri^/ml) ~ 1/40. Thus, our estimate 
for the induced value of #-term is: 

- A 5 -^^ logKK) * 3 • 1(T 13 . (24) 
8 47T 47r 

This value is just three orders of magnitude smaller than the current limit on the total 
#-term. Its contribution to the EDM of neutron could make sense only after specifying 
a mechanism of CP-strong puzzle solution. It is clear that the popular elimination of 
^-dependence due to axions makes the EDM of neutron unfeasible for (0). There are, 
however, some alternative solutions for CP-strong problem which assume the ^-relaxation 
at tree level only. This solution makes the most important the #-term contribution to the 
EDM of neutron. 

The last operator of interest is the Weinberg operator. Its distinguishing feature is 
the existence without hard gluon radiative corrections taken into account |J. However its 
numerical contribution to EDM in the four-family case is unlikely to exceed 10~ 33 e. cm. 
The additional smallness here is related to the strong suppression from the renormalization 
to hadronic scale. 



5 The limit of infinitely heavy rrih and m 
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We now check if the estimate (|22 ) is correct enough and prove the absence of exact cancel- 



lation or additional suppression at three loop level. 

The expression of interest is a three loop integral where quarks and SU(2) gauge bosons 
with comparable masses are involved. To simplify the problem, we impose on masses an 
artificial condition 

m L m l > m t > m w- (25) 
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which would allow to use an effective scale separation. This limit in general breaks the per- 
turbative unitarity and forces us to take into account further electroweak loops. However, 
it may be used to estimate the first coefficient of perturbative series. In our case it brings 
some important simplification to the problem of the loop calculation. The most important 
contributions are then associated with longitudinal parts of W-boson's Green functions. 



Indeed, if the limit fl25|) is held, all characteristic loop momenta could range between m t 
and heaviest masses. From that scale, according to our assumption, m w could be regarded 
as a small mass and l/m^ in the longitudinal parts of W-propagator makes its relative 
contribution to the effect to be enhanced in comparison with that of Feynman parts. Thus, 
it is clear that the quantity of interest, d s , will be proportional to the factor mf/m^ and 
we omit possible contributions of order and 1/mf. In other words, it means that in 

t'Hooft - Feynman gauge we take into account only diagrams with charged scalar Higgses 
and neglect those with W-bosons. The accuracy of the assumption for real massive 
parameters is presumably about 1/4 if we keep masses of heaviest quarks around 500 GeV 
and it is sufficient for our purposes. 

Two different possibilities of W-bosons attachment to the fermion line are depicted in 
Figs, la and lb. These skeletons should be dressed by one hard gluon loop and external 
soft gluon leg. Our calculation of the CEDM is based upon the external field technique 
proposed by J. Schwinger for QED and then extended on to the QCD case by Novikov, 



Shifman, Vainshtein and Zakharov (see for ex. the review |12[)- It can surely be applied 
to this problem because all loop momenta are much larger than the characteristic hadronic 
scale. The technique deals with the operator 

(x\Pp\y) = (x\iD,\y) = (i A + g s ^A^x))5\x - y) (26) 

where A c ^(x) is the gluonic field. Then the quark propagator taken in the background 
gluonic field reads as: 

(vac\Tq a (x)q b (y)\vac) = (x, a\i(P - m)' l \y,b) = (x,a\(P-m) , — d^' 6 )' 

P + Wsl 2(Gcr) — m z 

(27) 

where P = The field strength originates here as a result of commutation of two P's: 



[P„,P v ]=ig a G% / ^=ig a G tiI ,. (2* 



Using this technique, it is easy to demonstrate the vanishing of CEDM at two loops. 
This cancellation occurs before the last integration over a momentum of outer W-boson 
loop. The part of the amplitude of interest, Fig. la, comprises mass operator between two 
propagators antisymmetrized in masses: 

P - P 

— - r ^2 : - ( m i m i)i ( 29 ) 

P — mj P — mf 



The sum over i and j is performed according the prescription fllq). Here T denotes the 
mass operator taken in the background gluonic field. It possesses V - A gamma matrix 
structure and allows for the expansion in series of external field operators of increasing 
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dimension with some invariant functions depending on P 2 as coefficients. The explicit 
antisymmetrization in masses in (|29| ) leads to the following expression for this amplitude: 



K 2 - ™ 2 ) — — 2 [f, (P) 2 ] — — 2 (30) 

3 (P- m j)(P- m f) (P -m 2 )(P -mf) 

The commutator in this expression could be calculated for all operators entering to T. It 
could be shown that the result of commutation starts from operators with several powers 
of field strengths or field derivatives. Therefore, the two-loop amplitudes cannot induce 
CEDM (see Ref.0 for details). Clearly, at three loop accuracy this commutator must be 
a source of hard gluon field which compensates some extra dimensions. Another point of 
hard gluon attachment is not fixed and all other fermion lines should be expanded up to 
the first order in its field. The similar procedure could be performed over diagrams of the 
second type. 

It is clear that the result of integration may contain some power of \og(m 2 h /m 2 ) or 
log (m 2 /m 2 ). Naively we can regard this logarithm as a big parameter and calculate all 
diagrams in the "Leading Logarithm Approximation". To this approximation we believe 
that \og{m 2 Jm 2 ) <C ^ log 2 (m 2 /m 2 ) <C |, log 3 (m 2 /m 2 ). It is not true for real values of our 
masses. However to obtain an order of magnitude, we put all h log n equal to 1 in the final 
answer. The reason for assuming this artificial condition is to simplify the set of multi-loop 
calculation reducing it to the consequence of factorized integrations. 

It turns out that the logarithmic accuracy allows to find all relevant operators in T 
expansion and single out leading diagrams for the rest of amplitude. We start from the 
smallest distances determining the internal mass operator F for diagrams in Fig. la and 
one of mass operators in Fig. lb. These distances are associated with the propagation 
of g and h quarks respectively. The V - A gamma-matrix structure together with gauge 
invariance give three possible operators: 

r = — (c P 3 + cuvgsDpGr + c 2 g s ^{P u , G Q/3 }e^ Q/3 ) ^—-^ = (O + O x + 2 ), (31) 
The choice of these operators is not unique. We could use another basis; {P,P 2 }, for 

» 3 

example instead of P , etc. However, the expansion ([31]) appears to be the most convenient. 
The first term in this series, Oq operator, identically vanishes being substituted into fl30l) . 
For the second type of diagram this operator gives m 3 being applied to the s-quark wave 
function and therefore can be omitted as well. Two other operators in this series enter here 
with numerical coefficients c\ and c 2 whith no momentum or mass dependencies. They 
could be calculated from the expansion of mass operator in P at P 2 C m 2 . 

- if = /-4kg 1 - 75 , 1 P—P—P—* 1 -^ (32) 

2m^ J (2n)' l q 2 2 q — m q — m q — m q — m 2 
We used here the W-boson's Green function in the unitary gauge: 



: 9»» ~ QnQu/m 2 w q^q u 

2 2 — 2 2 V / 

q l — m l m mfcl 
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It is easy to see that the logarithmic divergence of (|32"D, converting to log(m 2 /p 2 ) after 
applying GIM mechanism, is associated with the Oq operator only The straightforward 
calculation leads to the following result: 

Ci _OL w 1 C2 &w 5 

We are left with two-loop expressions which could give a square of logarithm in the 
final answer. This means that the integration over a momentum flowing at the t quark line 
should be performed in last turn as we integrate from small distances to large ones. The 
corresponding diagrams and top quark propagation are shown in Figs. 2a and 2b. The 
blob represents operator T; dashed line here is the hard gluon propagator. It is clear that 
to logarithmic accuracy, the result of hard gluon loop integration may be presented as an 
effective vertex of W-boson with a fermion line changing the flavour from t to s. 

The position of external gluon field is not indicated in Figs. 2 and all propagators 
should be taken in the background gluonic field. At first look, the perturbative expansion 
breaks the main advantage of the calculation in the external field - we have to fix explicitly 
the gauge of SU(3) field. This problem could be resolved by dividing the four-potential A® 
into two parts: 

K = (AlU + a«, (35) 

where (A^) ext is the vacuum field while denotes the hard gluon field. It is a matter of 
convenience to choose an additional term fixing the gauge of a° in the following form (so 
called background gauge): 

-~WX) 2 > (36) 

where D e * l a a v = d^a^ + g s f abc {A b ^) ext a a v . It is easy to see that the gauge invariance for 
(Afyext field is still preserved. The hard gluon propagator in the external field in this gauge 
takes the form [ 12| : 



J zexp{ iq x}d*x(Ta«(x)a b u (y)) = (6 ab ^ - 2ig s f°* b A c x (y)^ + 2 -^r b G% + ...) exp{*g*}, 

q q q (37) 

where we have omitted the subscript "ext". This form of propagator is very useful for 
our problem and it allows us to calculate CEDM covariantly without fixing the gauge of 
external field. 

To reduce the number of diagrams the limit of large N c is used. In that limit combina- 
tions t a t a t b ; if abc t a t c ~ yt 6 are much larger than t a t b t a = -^t b . It means in particular that 
the external chromoelectric field can not be attached to the fermion line inside hard gluon 
loop as it shown at Fig. 3. The accuracy of this approximation is presumably -^2=1/9 and 
it is within the errors connected with previous assumptions. 

Let us first evaluate the contribution to this vertex from the operator 0\ ~ / y u D ll G fll/ 
which is a usual "penguin". By virtue of the equation of motion 

D,G% = -g s q^t a q (38) 

and Fiertz transformation this operator generates a close fermion loop as it is shown in Fig. 
4. The use of this equation is valid if the characteristic momenta inside this loop could be 
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regarded smaller than mh(m g ) which is satisfied to logarithmic accuracy. Both topologies of 
W-boson attachment generates a nonvanishing contribution to CEDM of s-quark. It turns 
out, however, that to logarithmic accuracy there is a cancellation between these two types 
of diagrams and therefore penguins operators cannot contribute to the effect at the level 
of log 2 {m\/m 2 ). It is possible to demonstrate that there is a nonvanishing contribution 
from these operators beyond this logarithmic approximation with a dependence of m 2 h jm 2 g 
parameter. It requires tideous true multi-loop calculations which is beyond the scope of 
our purpose. 

A similar cancellation between two topologies occurs when we substitute O2 operator 
and look for the effective vertex of W-boson with fermion. This operator, however, generates 
additional diagrams which should be taken into account as well. The first one represents a 
two loop mass operator inside W-boson loop (Fig. 5a); the second is the flavour changing 
mass operator inside gluonic loop (Fig. 5b). The crosses here indicate the chirality flips on 
the fermion line. 

The calculation of inner mass operator depicted in Fig. 5a in use of operator O2 is 
simple. To logarithmic accuracy, it is given by 

c 2 a s N c m 2 ^ 1-75 , QQ v 

M = ——, — rn t log(-[)g s {Ga)——, (39) 

7774 47T p A 2 

where p is the momentum flowing at the fermion line. It ranges between m t and the 
lightest mass of fourth generation quarks. The summation over flavours annihilates all 
other structures without the chirality flip. The only structure which does not vanish is the 
operator ([39|). Let us assume for the moment that m\ <C m 2 . Then this sum takes the 
form: 

77l 2 

{(g*),p) [ ^ iogg) - io g (-f )+ 

,m 




l(Ga),p]log(^), (40) 



(p 2 — m 2 )(p 2 — m\) m\ J ' p 



where we have neglected a further noncommutativity of momenta resulting from operators 
of dimension higher than that of CEDM. It is very natural that the logarithm is cut off 
at the lightest mass between m^ and m g . The last integration is also trivial and gives the 
CEDM operator of s-quark with a coefficient 

G F a s N c 1 c 2 m s m 2 1 m^m 2 ,) 
- I™{ V t s VtbV hb V hs )-^ — — - „ - log ( " )■ (41) 



v/2 4tt 16tt 2 m 2 w 2! 



mi 



The last graph given by Fig. 5b contains flavour changing mass operator. Its value 
at the incoming momenta k, which is much larger than m t , before the renormalization is 
given simply by: 

M(k) = kf(m - 75 ) = Sdjj* log(£)i^» (42) 

where A is the cut-off. The on-mass-shell renormalization prescription with respect to 
different masses to the left and to the right from the mass operator could be found in Refs. 
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0, 0. In our case this prescription looks trivial because both these masses are negligibly 
small in comparison with electroweak scale: 

M r (k) = k{f{k 2 ) - f(k 2 = o))(i - 75 ) = ~^i|r^^)^p. (43) 

The remaining integration over k ranging between m t and mh{m g ) is equivalent to that 
performed several steps earlier, in ([39]). The corresponding contribution to CEDM of s- 
quark differs from (|TJ]) by the factor -3. Combining these two numbers together and taking 
into account the value of the coefficient C2, we obtain the final answer for the CEDM of 
s-quark for very heavy fourth generation, to double logarithmic accuracy and in the limit 
of large N c : 

~ G a s a w 5N c m 2 1 2 ,m 2 (m 2 ) 

d s = -MV ts V tb V hb V hs )^rn s — — —yo g (— -^). (44) 

Substitution of N c = 3 and ^ log ( h m2 9 ) ~ 1 to this formula yields the estimation of 
CEDM close to 



6 Conclusions 

We have demonstrated a new interesting feature of the model with four generations of 
quarks incorporated into the same KM mechanism. The enhancement of neutron EDM in 
comparison with SM prediction comes from small distance effects which provide a regular 
factor of order m 2 / m 2 . The resulting value of EDM in this model, however, is just two orders 
of magnitude larger than corresponding SM value. The reason for that is in the numerical 
importance of large distance contribution in SM which is two orders of magnitude bigger 
than that coming from EDMs or CEDMs of quarks. 

The estimation of relevant operators in the effective Lagrangian performed in this work 
allows one to consider the influence of the fourth generation on other low-energy CP- 
violating observables such as T-odd form-factors of heavy nuclei. This problem deserves 
special consideration. 
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